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ABSTRACT 

Relativistic electrons moving into a highly tangled magnetic field emit jitter radia- 
tion. We present a detailed computation of the jitter radiation spectrum, including 
self-absorption, for electrons inside Weibel-like shock generated magnetic fields. We 
apply our results to the case of the prompt and afterglow emission of gamma-ray 
bursts. We show that jitter emission can reproduce most of the observed features with 
some important differences with respect to standard synchrotron, especially in the 
frequency range between the self-absorption and the peak frequency. We discuss the 
similarities and differences between jitter and synchrotron and discuss experiments 
that can disentangle the two mechanisms. 
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1 INTRODUCTION 

The external shock model (Meszaros & Rees 1997; Piran 
1999) has been very successful in the explanation of after- 
glow radiation in gamma-ray bursts. In this model the af- 
terglow photons are produced by converting some of the 
internal energy of the burst blastwave into radiation. The 
radiation mechanism is supposed to be synchrotron pro- 
duced by a population of relativistic electrons gyrating in a 
high intensity magnetic field. Synchrotron radiation is able 
to reproduce the spectra of some observed GRB afterglows 
(Wijers et al. 1999; Panaitescu & Kumar 2001; Panaitescu 
2005), their temporal decays, their duration and their small 
level of linear polarization (Covino et al. 1999; Lazzati et al. 
2004). 

How the magnetic field is produced is, however, a mat- 
ter of open debate. Compression of the interstellar magnetic 
field would generate a field that is not large enough to ex- 
plain the observed frequencies. On the other hand, a strong 
magnetic field linked to the central object powering the ex- 
plosion (e.g. the magnetic field of a neutron star) decays too 
rapidly with radius. In the standard external shock model, 
it is assumed that a quasi-equipartition magnetic field is 
generated by the shock and permeates the shocked inter- 
stellar medium (ISM) region without decaying. The frac- 
tion of energy given to the magnetic field is usually defined 
through B 2 = 8tyeb p where p is the energy density of the 
post-shock material and eg is a non-dimensional parame- 
ter called the magnetic field equipartition parameter. Syn- 
chrotron radiation is also supposed to generate the prompt 



emission photons. The synchrotron interpretation of prompt 
GRB spectra has however been unable to explain a sizable 
fraction of the spectra that display low-energy slopes steeper 
than v 1 ^ 3 , i.e., harder than what is allowed by synchrotron. 
In addition, models in which the electrons are accelerated 
impulsively suffer from the cooling problem, i.e., while the 
cooling time of the electrons is extremely short, no sign of an 
aging electron population is observed in the spectra (Ima- 
mura & Epstein 1987; Crider et al. 1997; Preece et al. 1998; 
Ghisellini et al. 2000). In this paper we do not attempt to 
solve the cooling problem, and implicitly assume that it can 
be solved by allowing for a slow or repeated acceleration of 
the emitting electrons. 

Particle in cell (PIC) simulations of collisionless shocks 
have demonstrated that the Weibel instability can produce 
magnetic fields at almost equipartition (e_g ~ 0.1, Silva et 
al. 2003; Fredriksen et al. 2004). Recent state-of-the-art PIC 
simulations demonstrate that the Weibel-generated mag- 
netic fields survive for at least a few hundred skin depths 
(Spitkovsky 2005; Chang, et al. 2007; Spitkovsky 2007), 
in agreement with predictions by Medvedev, et al. (2005). 
Whether these fields can avoid dissipation and survive be- 
yond the few hundred plasma skin depths remains an open 
question. What appears certain is that the field created by 
the Weibel instability has a very short coherence length, so 
that classical synchrotron formulae cannot be generally ap- 
plied. 

Radiation produced by relativistic particles in small- 
scale magnetic fields has been a subject of study for several 
decades. Historically, it begins with two seminal papers by 
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Landau & Pomeranchuk (1953) and Migdal (1954), who pre- 
dicted the suppression of the spectral power of low-energy 
harmonics — the Landau- Pomeranchuk-Migdal effect. In- 
terestingly, the spectral properties of synchrotron radiation 
from either homogeneous fields or from random fields co- 
herent on scales larger than the particle Larmor radius is 
rather universal, as it can be straightforwardly applied in 
a wide variety of magnetic field configurations. In contrast 
to synchrotron radiation, the radiation emitted from a mag- 
netic field (not necessarily random) whose gradient scale is 
smaller than the particle gyro-orbit, is determined by the 
particular field configuration. Therefore, one has to be very 
specific about the field structure in the system for which the 
radiation spectrum is computed. 

A number of theories of radiation by relativistic elec- 
trons in small-scale magnetic fields in various regimes have 
been developed, including radiation in wiggler/undulator 
fields and free-electron lasers (presently used in sources of 
X-ray and extreme UV light) (Kincaid 1977; Joshi et al. 
1987; Fedele et al. 1990; Williams et al. 1993; Attwood et 
al. 1993, Attwood 2000), radiation emitted by cosmic rays 
in the interstellar medium plasma turbulence if it extends 
to sufficiently small scales, beyond the Larmor radius of 
thermal electrons (Kaplan & Tsytovich 1969; Bel'kov et 
al 1980; Ginzburg & Tsytovich 1984; Toptygin 1985; Top- 
tygin & Fleishman 1987ab); transition radiation in inho- 
mogeneous plasmas (for a good review, see monograph by 
Ginzburg & Tsytovich 1984); synchro- Compton (or nonlin- 
ear inverse Compton) radiation produced by electrons in the 
non- stationary field of a strong electromagnetic wave (Gunn 
& Ostriker 1971; Rees 1971ab; Blandford 1972); betatron ra- 
diation emitted by a particle propagating in a magnetic field 
of an ion current channel (somewhat analogous to a Weibel 
current filament; Wang, et al 2002). 

Medvedev (2000) considered radiation from Weibel- 
generated fields, generalizing the wiggler/undulator radia- 
tion to the case of randomly distributed current filaments 
and associated fields as a model of a GRB shock. The spec- 
tral properties of radiation have been analytically derived, 
using the perturbation theory (Landau & Lifschitz 1971), 
for a single electron and for a power-law electron energy 
distribution, in the limit of a magnetic field with a very 
small coherence length. The theory of jitter radiation ap- 
plied to the prompt emission of GRBs was shown to be able 
to reproduce some characteristic features, such as the steep 
low-energy spectrum and the sharpness of the spectral break 
(Medvedev 2000). The theory has been further generalized 
to the 3D case by Fleishman (2006b) who calls it "Diffusive 
Synchrotron Radiation" (DSR) and by Medvedev (2006). 
Jitter radiation and diffusive synchrotron radiation describe 
the same phenomenon: radiation from particles traveling 
through small scale, random magnetic fields coherent on 
the scale of a plasma skin depth. The primary difference 
in the two approaches is that Medvedev uses the pertur- 
bative approach introduced by Landau & Lifschitz (1971) 
whereas Fleishman applied both the perturbative treatment 
and a more general non-perturbative scheme which allows 
for the effects of scattering to appear in the spectra as the 
coherence length of the field increase. Fleishman (2006b) 
and, independently, Medvedev (2006) obtained the single 
electron spectrum to be F v oc u° in the perturbative ap- 
proach, noting the possibility that steeper spectra (up to 



v ct v 1 ) could be obtained if the magnetic field spectrum 
could be factorized. Previous results were obtained by solv- 
ing the non-perturbative equation describing single electrons 
emitting radiation due to both large scale ordered fields and 
small scale random fields (Toptygin 1985; Toptygin & Fleish- 
man 1987a). The anisotropy of jitter radiation has also been 
suggested as a possible explanation of the rapid spectral 
variability of prompt GRB emission (Medvedev 2006). Re- 
cently, the jitter radiation model has been extended to the 
self-absorption regime (Medvedev, et al. 2007) and applied 
to GRB afterglows. Approximated analytical spectra and 
light curves were computed. 

In this paper we use the theory of jitter radiation as 
it is more amenable to numerical and analytical calcula- 
tions (following Medvedev 2006) to improve upon the an- 
alytical computations of Medvedev et al. 2007. We review 
and refine the jitter theory including self-absorption, numer- 
ically implement it in a radiation code, and apply it for the 
optical-UV range for the prompt emission of GRBs. We also 
compare properties and spectral characteristics of jitter and 
synchrotron radiation models in the afterglow regime. We 
show that jitter radiation is a viable framework for the inter- 
pretation of prompt and afterglow observations. We discuss 
the observational implications of this alternative radiation 
mechanism and discuss some observational tests to compare 
synchrotron and jitter radiation in GRBs. 

The paper is organized as follows: in § 2 we compute 
the emissivity of jitter radiation and the absorption coeffi- 
cient, in § 3 we discuss the range over which our results are 
valid. In § 4 we apply our result to some cases relevant for 
GRB prompt and afterglow emission. We conclude in § 5 by 
discussing some implications for GRB observations. 



2 SPECTRUM OF JITTER RADIATION 

This section lays out the analytical framework we use to 
compute the specific intensity 1^ of jitter radiation produced 
by Weibel generated magnetic fields. In the first part we de- 
tail how we mathematically describe the magnetic field pro- 
duced by the Weibel instability and used in the calculations. 
In the second part we discuss the mathematical formulation 
used to calculate the emission coefficient , which describes 
the angle- averaged power per unit frequency added to the 
radiation field by electrons emitting jitter radiation. Finally, 
the third part describes how to obtain the absorption coef- 
ficient a^, which describes how jitter radiation self attenu- 
ates. 

As a first step we derive the absorption and emission 
coefficients in the comoving frame. Although we calculate 
the comoving quantities, it is straightforward to transform 
to the observer frame by recalling that a^co and are 
both Lorentz invariants with lo — 0/7;, (1 + \ cos 6') where 
7t and v refer to the bulk velocity of the object and not 
to the Lorentz factors of the individual emitting electrons 
and cos Q' is defined by the angle between the photon path 
and the co-moving volume velocity direction. This algebraic 
transformation allows us to recast the equation for the ob- 
served intensity as a function of comoving quantities. Pro- 
ceeding in this manner greatly simplifies the analytics and 
allows us to solve for the observed specific intensity by using 
numerically computed, comoving quantities. 
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2.1 Magnetic Fields In Gamma Ray Bursts 

The standard model for emission from both the prompt 
burst in GRBs and the afterglow relies mainly on the 
assumption that the radiative mechanism is synchrotron 
(Zhang & Meszaros 2004; Granot et al. 1999ab). While the 
assumption of synchrotron is reasonable for some of the ob- 
served prompt bursts due to strong magnetic field progen- 
itors, it is unlikely that it is always dominant in afterglow 
emission. Parameterizing the strength of the magnetic field 
by €b, which is the ratio of the magnetic field energy den- 
sity to the thermal energy density, leads to values of es 
(from observations and simulations) ranging from 10 _1 to 
10~ 5 (Medvedev & Loeb 1999; Panaitescu & Kumar 2001; 
Panaitcscu 2005). 

The field may be directly connected to the compact 
object powering the GRB outflow. In that case, assum- 
ing a strong field of ~ 10 16 G on the surface of the com- 
pact object and considering volume expansion in an outflow 
(B oc V~ 2/3 ), one obtains B ~ 10~ 2 G (e s ~ 10~ 7 ) in the 
afterglow region (R ~ 10 16 cm). Such a field would be too 
small to reproduce the observations and would decay too 
steeply with time. Alternatively, the magnetic field could 
be generated through shock compression of a pre-existing 
interstellar field. In that case, Bl oc 'yBisM is too weak 
to generate sufficiently strong fields at the afterglow shock 
front and results in values of es ~ 10~ n (Medvedev & Loeb 
1999). With no new mechanism available to produce fields 
it is unlikely that synchrotron can always be the mechanism 
by which afterglow radiation is produced. 

One very promising candidate for the origin of the mag- 
netic field which results in afterglow emission comes from a 
relativistic version of the well known two-stream Weibel in- 
stability in a plasma. In the rest frame of a relativistic shock, 
the Weibel instability amplifies any existing magnetic field 
perturbations by growing current filaments out of instream- 
ing electrons. The amplification of the current filaments is 
fed by the kinetic energy of the inflowing material (the ISM 
in the case of GRBs) and does not saturate until all of the 
energy in the particle distribution function anisotropy is con- 
verted into magnetic field energy. Medvedev & Loeb (1999) 
have shown that the field is generated in the plane of the 
shock and has a coherence length of order the relativistic 
skin depth of the shock. Further numerical studies (Noshed 
et al. 2003; Silva et al. 2003; Fredriksen et al. 2004) have 
confirmed the existence and growth of this instability. For 
a detailed discussion of the Weibel instability see Medvedev 
& Loeb (1999). The major factor that makes magnetic fields 
generated via the Weibel instability unique is that the cor- 
relation length of the fields is less than a Larmor radius and, 
as a result, the radiation generated by electrons in Weibel 
fields is very different from the one generated by uniform 
large scale fields which result in synchrotron radiation. Ad- 
ditionally, the field is generated in the plane of the shock 
and is then transported downstream as the shock propagates 
into the ISM. This continual generation of field is an attrac- 
tive solution to the problem of how field is either generated 
or carried far enough downstream to produce the radiation 
observed in afterglows. 

In this paper we adopt the model introduced by Fleish- 
man (2006b) and Medvedev (2006) to describe the mag- 
netic field in Fourier space. We emphasize that such a 



parametrization of the field by no means incorporates all 
the detail and phenomenology of the fields observed in PIC 
simulations of collisionless shocks. However, it has the great 
advantage to be analytically tractable and relatively simple. 
Direct computations of radiation spectra can be coupled to 
simulations (e.g. Hededal 2005) but they are time consum- 
ing and, albeit with a very high accuracy, provide spectral 
calculations valid only for that particular simulation. We 
further note that we have explicitly excluded the presence 
of electric fields in our work. Hededal (2005) found that the 
electric fields are generated in PIC simulations and have 
~ 10% of the energy associated to the magnetic fields. He 
also concluded that the primary effect of electric fields was 
to flatten the low frequency slope. We expect the electric 
fields to play a secondary role in the generation of the emer- 
gent spectrum at the wavelengths we consider and as such 
we have simplified the analytics and numerics by excluding 
it. 

During the initial stage of the Weibel instability, the 
field is generated in the plane of the shock, with k± ~ 
skin depth and k\\ = 0. At this point, the parallel and 
transverse spectra are decoupled. Soon after saturation, sec- 
ondary instabilities develop in the parallel direction thus 
making k« > 0. Since the instability develops in the parallel 
direction only, there is no coupling to the transverse dynam- 
ics either, at this time, so the spectra are independent. Some 
coupling of the transverse and parallel dynamics begins to 
develop at late times, when the filaments start to twist, in- 
teract with each other and merge. At this point, we shall 
expect that our assumption of the separability of the field 
spectrum becomes an approximation. This is important es- 
pecially in the afterglow phase, since it has been shown that 
afterglow radiation has to be produced in the whole blast 
wave and not only behind the shock (Rossi & Rees 2003). 
However, as long as the current filaments are not completely 
intertwined to form a statistically homogeneous state, one 
shall expect some degree of statistical independence of the 
field spectra to remain. In our description, the correlation 
tensor for the field is given by 

if a/3 (k) = 4nC(8 al 3 - n a n0)f z (k\\)f xv (k±). (1) 

We assume that the field in the plane of the shock and 
the field in the direction of propagation can be factored from 
each other where f xy (k± = (fc 2 + fc 2 ) 1 / 2 ) describes the field 
in the plane of the shock, f z {k\\ = k z ) describes the field 
in the direction of the shocks propagation, and C is a nor- 
malization constant. The functions used to describe the field 
distribution are taken from Medvedev (2006) in order to fa- 
cilitate a direct comparison with the results therein and are 
given by 

f '^ ) = J^T^' (2) 

and 

k 2ct2 

fxy{k±) = k + k\ (3) 

where cti,2, /?i,2, and k ± y are parameters used to fit the 
spectrum to numerical results. In this paper we have chosen 
ai = Q2, Pi — @2, and k± — k\\ = ks of the local field. 
In general k is a parameter determined by local quantities 
(Medvedev et al. 2005) but is treated as a constant in this 
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paper as it does not significantly alter the shape of the spec- 
trum in the regimes under consideration. 

Finally, C is fixed by the requirement that 
/ K aa (k)dk = (B 2 ) where (B 2 ) is the mean square 
value of the local magnetic field . This convention for C 
results in a normalization of 



C 



(B 2 ) 



ffz(k\\)k r f r (k r )dk^dk r 



(4) 



where we have switched from Cartesian to cylindrical coordi- 
nates to simplify the integration. Unlike previous works, we 
have chosen to normalize the correlation tensor by the mean 
square value of the magnetic field in order to obtain abso- 
lute values for the luminosity of the afterglow and prompt 



2.2 Emission - j u 

We now calculate the emissivity of an ensemble of isotrop- 
ically distributed electrons in the jitter regime. The emis- 
sivity of a power-law distributed electrons, N oc 7 _p with a 
sharp low-energy cutoff 7 < 7min has been calculated for the 
simple one-dimensional jitter model in the original paper by 
Medvedev (2000). In two separate, fully three-dimensional 
treatments of jitter radiation (Fleishman 2006b; Medvedev 
2006) the single electron spectral power has been calculated. 
The ensemble emissivity is computed as a convolution of the 
single electron spectral power with the electron distribution. 
We present here, for the first time, the total radiation emit- 
ted by a distribution of electrons using the coefficients which 
return the true spectrum as a function of local conditions. 
Other than the powers on the correlation tensor describing 
the magnetic field (to which the radiation spectrum is rela- 
tively insensitive) every attempt has been made to keep the 
number of free parameters to a minimum. 

The formula used to describe the angle averaged ra- 
diation emitted by a single, relativistic particle traveling 
through small scale magnetic fields is given (neglecting 
plasma dispersion which change the low frequency results 
to oc uj 2 ) by (Landau & Lifshitz, 1971, section 77 p. 215) 



dW _ e 2 u) 

diU 2-7TC 3 



1 - 



ui 7 



2lj' 2 j 4 



du', (5) 



where 7 is the Lorentz factor of the particle. Due to ap- 
proximations used in deriving equation ([5}, this equation is 
only valid when the ratio of a particle angular deflection due 
to magnetic field fluctuations (a) to its relativistic beaming 
angle (A8 ~ I/7) is much less than one. This approxima- 
tion is expounded upon later in the paper. In keeping with 
Medvedev (2000) we define this ratio as 



kspe 



a 

~A8 



(6) 



where p e is the Larmor radius of an electron. In practice the 
value of 8 calculated, for our parameters, is ^ 1. The value 
of 8 enters into the integration limits in Eqs. llOl and llll below. 
and setting it equal to zero there has a small effect. The effect 
of this assumption can be seen in Figure 4 of Medvedev 
(2000). Intermediate values of 8 produce an upturn in the 
spectrum just below the peak that disappears when 8 0.3. 
The shift in the peak frequency and the normalization, due 
to the change in 8, are not lost in our computation. 



The term |w^/| 2 in equation ([5)l is the square of the 
Fourier transform of the acceleration field due to the Lorentz 
forces. Here we replace it with a volume averaged (lw^/1 2 ) 
by assuming a statistically homogeneous turbulence. The 
derivation of this term is left for the appendix and has been 
included for the sake of completeness (the interested reader 
is referred to Fleishman 2006b for the full details) and we 
state here the result assuming the fields described above 
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where 1(0') is given by 



= |/ ! ft|)/.«(*iW«' + k'v) d 3 k, 



(7) 



(8) 



where O' is the angle between the particle velocity and the 
observer in the comoving frame, C is given by equation Q, 
T is the period for an electron traveling in Weibel fields, 
and fz(k\\) and fxy(k±) are given by equations (J2j) and Q, 
respectively. To evaluate the integral in equation (Jsj) it is 
necessary to specify both the limits of integration and the 
geometry of the problem. The limits come from the nature 
of the magnetic fields generated by the Weibel instability. 
The Fourier component associated with the fastest growing 
mode in the Weibel instability is 



47shockW pe 



V2, 



(9) 



The factor of 4y s hock comes from the shock compres- 
sion and uip e — 4ne nga:t is the plasma frequency of the 
pre-shocked material (which, in the case of afterglows, cor- 
responds to the ISM or wind density profile). While modes 
are initially compressed only perpendicularly to the shock 
plane, we make the assumption that the spatial scales are 
mixed by turbulence and set the inverse length scales, fcy, 

kj_ = kweibel- 

The geometry of the system is defined as follows: The 
shock is propagating in the z direction and lies in the x-y 
plane. A particle in the shock has a velocity vector given by 
k = xfcsinO' + zfccosO', this gives us k • v = fc^wsinO' + 
kzVeosO' . Medvedev (2006) presented three separate forms 
for (lw^/1 2 ), corresponding to a shock viewed at degrees, 
at tt/2, and in between these extremes. We chose to evaluate 
equation (J5J somewhat differently by using the properties of 
the delta function to integrate the correlation function in two 
distinct ways. By doing this we can match the asymptotic 
forms with two functions as opposed to three and do not 
suffer from numerical errors. The two forms are necessary 
to avoid the introduction of infinities as O approaches and 
7r/2 degrees. The two forms we use are 

M9') = J |^ei/.(^7 + fetanG') JW(M^,(10) 
and 



77 fz(kz)fxy ( . + 



\v sin O 



k z 



sin 6' tan 6' 



dk z , (11) 



where the bar in f X y denotes an integration over k y . In 
Fourier space, the modes that generate the field lie within a 
spherical annulus with radius k± and ku G [8k Weibe i, k Weibe i] 
and this sets the limits for the integrals in equations (|10|l 
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and Our magnetic field model assumes that fluctu- 

ations in the magnetic field are not generated on physi- 
cal scales larger than ~ (Skweibei)' 1 and therefore these 
scales are excluded from the integration region. In practice, 
these large scales contribute little to the overall spectrum 
for S <C 1 and this is seen by the lack of any change when 
we replace Skweibei by zero in our code. 

Choosing v ~ c and switching between forms for 0' ~ 
j we are able to reproduce the results of Medvedev (2006) . 
Combining equations ©, ©, ©, ©, ©, JTOj, and (JTTJ) 
and multiplying the result by ^ gives us 



y ' T duj 



(12) 



the total energy per unit frequency and unit time, emitted 
by a single electron in a magnetic field with the configu- 
ration specified by Eq. [T] Unlike the case of synchrotron, 
where it is possible to precisely define an orbital period, the 
random, small scale nature of Weibel turbulence requires a 
more arbitrary choice for T. We follow Medvedev (2006) and 
Fleishman (2006b) which includes the value of the period in 
the magnetic field derivation (see equation @). By using 
this method, T simply disappears from the final form for 
P{u). 

The final step required to calculate the emission from a 
population of electrons is to choose the form for the distri- 
bution. It is standard in GRBs and other sources of non- 
thermal radiation to assume that the electrons follow a 
power law distribution 



n (7) = Kj p 



where 7 G [7mm, 00], K = (p — l)n' 7 ; 



1 ~p— 1 



and 



2 e e e' 
1 n'm e c 2 



(13) 



(14) 



The primed quantities n' and e' refer to the comoving num- 
ber and energy densities and e e is the fraction of the thermal 
energy in the electrons. 

Combining equations (|13[) and (|12|l and integrating 
from 7 m i„ to 00 yields the final result 



Ptot(w) 



n ' { / y)P{io)d r y. 



(15) 



It is this equation that we have numerically solved for 
several values of 0'. In order to compute j u we make several 
assumptions. First, we assume that there is no scattering. 
Second we assume an isotropic electron distribution. Third, 
we assume that for a given point on a surface, whether it 
is a plane parallel slab or a spherically expanding shell, we 
only see the radiation due to the comoving angle aligned to 
the local bulk velocity. This third assumption is physically 
equivalent to assuming we see only 'flashlights' from differ- 
ent emitting regions. Finally we assume that the medium 
is moving ultra-relativistically and is beamed into a narrow 
cone of angular width I/7. Using these assumptions we make 
the simplifying approximation that — P to t(u)S(£l — Q(r)). 
Bearing in mind that we have so far calculated all quan- 
tities in the comoving frame, this form for the emissivity 
introduces an additional factor of 7 2 (1 + /3cos0') 2 when 
transforming to the observer's frame due to the transform 
of the solid angle (Granot, Piran, & Sari 1999a; Rybicki & 
Lightman 1979). 



In order to check our results we have performed the 
following tests. First, it can be shown analytically (for 
/3i,2 = 0) that combining the 0' = form of (Iw^l 2 ) with 
P = (2e 2 7 4 /3c 3 )(™i + 7 2 «; 2 ) (equation 4.92 of Rybicki & 
Lightman 1979) Eq. PT21 reduces to 



dW/dt = (2/3)r 2 c7 2 (P 2 ) . 



(16) 



which is equivalent to the case of synchrotron radiation. We 
have confirmed that our procedure numerically returns the 
same value by integrating equation JJj). Setting /3i,2 7^ in 
our code does not change this result. Second, we replaced 
P(uj) in equation (I15|l with the the form for synchrotron ra- 
diation and verified that our technique yielded results which 
were in agreement with the analytic form for Ptot, Synch (u) 
given by equation 6.36 in Rybicki & Lightman (1979). 



2.3 Absorption - 

Given the single particle emissivity, P(v), and the parti- 
cle energy distribution, the calculation of the jitter self- 
absorption is straightforward. We have analytically calcu- 
lated it (Medvedev, et al. 2007) in the regimes when the self- 
absorption frequency is above and below the jitter peak fre- 
quency. For numerical implementation, the self-absorption 
coefficient can be derived using equation 6.50 from Ry- 
bicki & Lightman (1979) as follows: 



Sixv 2 



y ' ' dE E 2 



Using P(y,E) = 2tvP{u), E = ~fm e c 2 , N(E) - 



(17) 



^n( 7 ), 



and oj = 2nv where n( 7 ) and p(u)are given by equa- 
tions (|12[l and (|13l) . quickly yields 



(p + 2)7T 2 fc 



7 



P(u)dy. 



(18) 



Effectively, once the numerical problem of solving equa- 
tion (|15|) has been solved, the only modification necessary 
to obtain the absorption coefficient is a change in the power 
on gamma and a change in the normalization constant. In 
using this form for the absorption coefficient we have as- 
sumed that the electron distribution is locally isotropic, that 
no scattering occurs, and that the local absorption is domi- 
nated by electrons moving with only one angle with respect 
to the observers line of sight. For a further discussion of the 
absorption coefficient see Medvedev et al. (2007) . 



3 LIMITATIONS OF VALIDITY 

Equation ([5]) is derived using the perturbation theory (Lan- 
dau & Lifshitz 1971) and, consequently, has a limited range 
of applicability. Jitter radiation explicitly assumes that the 
coherence length of the emission is small enough that the ap- 
proximation of rectilinear motion is valid (Medvedev 2000 & 
Fleishman 2006b). The coherence length A c is proportional 
to 1/uj and clearly increases as the frequency decreases. 
Fleishman (2006) describes that, at lower frequencies, a par- 
ticle trajectory random walks due to repeated scatterings 
off magnetic field inhomogeneities. This scattering results 
in a modification of the emitted radiation in the region of 
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the spectrum just before where the plasma dispersion dom- 
inates. The non-perturbative version of the theory of radi- 
ation emitted in small scale field (Toptygin & Fleishman 
1987, Fleishman 2006a) accounts for this effect and finds a 
P(oj) oc w 5 regime in this region. At even lower frequencies, 
u><ujpi, the spectrum follows P(uj) oc uj 2 due to plasma dis- 
persion. This region is at photon energies of a fraction of 
an eV or less in prompt GRBs and at frequencies below 
~ 1 MHz in afterglows, assuming standard GRB parame- 
ters. 

The perturbative form of the theory fully accounts for 
the spectrum above and below (down to the plasma disper- 
sion regime) the P(ui) oc lu^ region. For the purposes of 
this work it can be shown that the region in frequency (ap- 
proximately) given by u < eBOJjm (where uij m is the peak 
frequency of jitter radiation for a given minimum Lorentz 
factor) is affected by scattering and therefore the spectrum 
obtained with the methods used in this paper is only ap- 
proximate. This limitation does not pose a problem for this 
work. It can be shown that S « (Medvedev et al. 

2007) which requires we assume e_e < 10~ 3 and, more re- 
alistically, in the range of 10" 6 - 10" 4 . Even for a value of 
es = 10~ 2 the cutoff frequency esiOjm corresponds to an 
energy at the edge of the BATSE cutoff window of 20 keV. 

The only case in which the limitations of the pertur- 
bation approach may be relevant is the modeling of the 
afterglow spectrum. Even in such a case, however, the re- 
gion in question generally lies entirely within or close to the 
optically thick portion of the spectrum, in which case the 
spectrum simplifies to the source function. In cases where 
the spectrum generated by our perturbative approach may 
be incorrect we will note it in the individual results. 



4 RESULTS 

We here consider a single internal shock within the prompt 
phase of a GRB. The material has a pre-shock density of 
9.25 x 10~ 12 g/cm 3 , a shock Lorentz factor of ji n t — 2, 
e e = 0.2, and a magnetic equipartition fractiorQ e_g = 10~ 5 . 
The spectral index of the electrons is p — 2.5. The thickness 
of the shell is 3 x 10 6 cm in the observer frame and the 
material is moving with a bulk Lorentz factor of 100. We 
consider a shell radiating at R = 3.2 x 10 14 cm. This setup 
corresponds to an isotropic energy of 5 x 10 52 ergs in the 
shell. 

These values are chosen to give a peak frequency of 
about 200 keV in the observer frame. The parameters de- 
scribing the magnetic field distribution do not affect the 
spectral characteristics as long as ai,2 > 0.5 and /3i,2 > 
oti,2 + p/2. To be consistent with Medvedev (2006) we set 
01,2 = 2. and f3\$ = 10. It is important to keep in mind that 
the emission angles correspond to co-moving angles. The re- 
sults plotted are the observed quantities as a function of 
co-moving emission angle. 

Figure [T] shows the emissivity vs. frequency for different 
emission angles (8'). For 8' = the emission has a slope 
q„ = 1 at low frequencies and then turns over to a slope 

1 We note that there is no experimental value for the 
value of eg during the prompt phase. 



a„ = — (p — l)/2 = —0.75, set by the electron power in- 
dex. As the angle increases, the low frequency slope near 
the peak decreases, and at much lower frequencies the slope 
flattens to a v — 0. Beyond an angle of 30°, there is no longer 
a peak in the spectrum, just a transition from a constant 
(a„ = 0) to OLv = —0.75. The value of the low frequency 
constant changes with angle, peaking at 45°. At high fre- 
quencies, the slope is always the same but the emissivity at 
a given frequency decreases with angle. These results are in 
full agreement with the single electron spectra of Medvedev 
(2006). Figure [2] shows the absorption coefficient a for the 
same setup. The absorption coefficient is proportional to 
P(uj) x lu~ 2 below the peak frequency and to P(lj) x lj~ 2 ' 5 
above the peak frequency. 

Below frequencies of ~ 3 x 10 -3 eV our perturbative 
approach is no longer valid. A non-perturbative approach 
would be needed to properly calculate radiation at lower 
frequencies (Fleishman 2006b). Therefore, frequencies in this 
regime are not considered in this paper. 

With P(uj) and a we can solve the ID radiative trans- 
fer equation to find the specific intensity of the shell at dif- 
ferent angles. Assuming the observable portion of the shell 
is spherical and that the amount of time any piece of the 
shell is radiating is short compared to the time it takes the 
shell to become visible (i.e., we assume the shell is infinitely 
thin), the specific intensity at a given angle is proportional 
to the total emission at the time that angle comes into view. 
Figure [3] shows the integrated specific intensity at different 
angles for jitter radiation. Figure [4] shows synchrotron radi- 
ation for identical conditions and angles. 

The synchrotron radiation spectrum shows no change 
in the shape of the spectrum. The peak frequency decreases 
by about a factor of 2 and the power decreases by a factor 
of ~ 16 between 0° and 90° due to relativistic effects. Jitter 
radiation, on the other hand, shows a large change in spec- 
tral index and total emission below the peak frequency and 
a decrease in emission of about a factor of ~ 27.5 with angle 
above the peak frequency. Note that the source function for 
jitter and synchrotron are the same for identical conditions. 
Self absorption becomes important at a few 10s of eV in this 
example. 

The time at which different angles come into view is 
found from the geometric time delay to be (1 — cos#) x r/c 
where 8 is the emission angle in the observer frame, r is 
the emission radius (10 14 cm in this example) and c is the 
speed of light. Figure [5] shows intrinsic brightness vs. time 
at 10, 100, and 1000 keV and 1 eV (124 nm). At 1000 keV, 
above the initial peak frequency, the emission decreases with 
time by a factor of about 46 between t — and t — 2.1. At 
100 keV, just below the initial peak frequency, the emission 
increases slightly for 0.1 s and then decreases. At 10 keV the 
emission increases by a factor of 7 in 0.25 s, corresponding 
to 40°, and then decreases. The initial amount of increase 
is larger at lower energies. At 1 eV, there is an initial sharp 
increase. Very early on, however, the emission becomes op- 
tically thick and begins a shallow decrease (by a factor of 
only ~ 7 with time). This behavior has important implica- 
tions for the interpretation of spectral lags and the optical 
emission during the prompt phase. As is clear from Fig. [5] 
the optical and high energy bands can show very little cor- 
relation, even if they come from the same electrons and the 
same emission mechanism. Lack of correlation was used as 
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an argument in favor of a reverse shock origin of the op- 
tical flash in GRB 990123 (Akerlof et al. 1999, Sari et al. 
1999). The flattening of the spectra in the X-ray band could 
also explain the soft X-ray excesses detected in some events 
(Preece at al. 1996; Strohmayer et al. 1998). 

Figure [6] shows the ratio of the total number of photons 
between 110 and 320 keV to the number of photons between 
55 and 110 keV vs. time. This ratio is equivalent to the count 
ratio of BATSE channels 3:2 (HR 3 / 2 ). This plot shows a hard 
to soft evolution for the first 0.4 s (out to a 45° angle) and 
then a fairly flat hardness ratio after this. 

Let us now consider the external shock phase. For an 
example comparison of the afterglow spectrum for jitter 
and synchrotron radiation we examine an isotropic explosion 
with an energy of 10 53 ergs expanding into a medium with 
a density of 1 particle/cm 3 at the observed time t — 1000 s. 
The Blandford-McKee (Blandford & McKee, 1976) solution 
for a relativistic fireball is used to determine the density 
(~ 100 cm -3 ), internal energy (~ 4 erg cm _3 )and Lorentz 
factor (F ~ 25) of the material at the leading edge of the 
shock. A thin shock model is used to calculate the emis- 
sion and the optical depth. A magnetic field equipartition 
parameter e_B = 10~ J and an electron equipartition param- 
eter e e = 0.01 are assumed everywhere, and a shock thick- 
ness of ~ 2.4 x 10 14 cm everywhere is assumed. Fig- 
ure [7] shows a comparison between the emission from jitter 
and synchrotron radiation under these conditions. Below the 
peak frequency, the synchrotron emission increases with fre- 
quency as v 1 / 3 , while the jitter spectrum is nearly flat. This 
is due to the fact that large angles dominate the afterglow 
spectrum (see also Medvedev et al. 2007). The synchrotron 
and jitter spectra become optically thick at around 10~ eV 
(24 GHz) and 10" 3 eV (240 GHz), respectively, in this exam- 
ple, but this is dependent on the thickness of the afterglow 
material. 



5 SUMMARY AND DISCUSSION 

We have computed jitter radiation spectra for a non-thermal 
population of relativistic electrons radiating in a highly non- 
uniform magnetic field, including self absorption. Our re- 
sults are in agreement, when relevant, with previous com- 
putations (Fleishman 2006b; Medvedev 2006; Medvedev, et 
al. 2007). Due to the complexity of the problem, several as- 
sumption and/or simplification were made in order to obtain 
a semi-analytical result for the spectra. We here review the 
assumptions and simplifications made. 

• The most important simplification we adopt is that of 
assuming that the field can be factorized in the parallel and 
perpendicular directions with respect to the shock propaga- 
tion velocity. This assumption is fundamental in our results 
as it allows for the dependence of the low-frequency spectral 
slope on the viewing angle. It is likely to be a good assump- 
tion immediately behind the shock and become progressively 
worse as the ion current filaments get twisted and/or merge, 
mixing the parallel and perpendicular components of the 
field. 

• The strength of the magnetic field is assumed to be uni- 
form throughout the shocked region and the average mag- 
netic field to be rather small (containing only a fraction 



~ 10~ 5 of the internal energy of the blast wave). Such as- 
sumption is necessary to simplify the computation and to 
make sure the radiation is produced in the jitter regime 
rather than in a transition region. PIC simulations (Silva 
et al. 2003; Fredrikssen et al. 2004; Cheng et al. 2007) find a 
stronger magnetic field right behind the shock and a moder- 
ate decline behind it. However, PIC simulations can follow 
the shocked material only several hundred skin depths be- 
hind the shock, while afterglow radiation has to be produced 
by the whole blast wave, in order to avoid a strong inverse 
Compton component that is not observed (Rossi & Rees 

2003) . 

• The peak of the magnetic spectrum can be around the 
ion skin length and this seems to invalidate the jitter condi- 
tion for leptons. However, this does not invalidate the use of 
jitter radiation, because the deflection angle depends on the 
product of the field correlation length and the field strength. 
We took this into account: the jitter break will be above the 
synchrotron frequency in the field of the same strength if 
es 10 -3 . Fits to observational data of GRBs favor such 
low values (Panaitescu & Kumar 2001; Table 3 of Panaitescu 
2005). If the field strength is larger, the jitter regime trans- 
forms into the wiggler regime (Attwood 2000) and, for even 
stronger fields and larger correlation lengths - into the syn- 
chrotron regime. We should also stress here that in reality 
the magnetic field spectrum spans at least several decades 
in fc-space, with the highest k being even higher than u> pe /c, 
even for the electron-ion plasmas (Wiersma & Achterberg 

2004) . This is also confirmed in PIC simulations. Radiation 
from these fields will surely be in the jitter regime; the low- 
k part of the field spectrum will produce synchrotron-like 
radiation, as is discussed in (Medvedev 2000). 

• Our computations are performed in the perturbative 
approach and cannot be applied in some regimes. Especially 
relevant is the low frequency range where the photons scat- 
ters off the magnetic field inhomogeneities resulting in a 
F(u) oc u 1/2 regime (Fleishman 2006b). 

• Simulation of magnetic field generation in ion-electron 
collisionless shocks find evidence of electric fields on top of 
magnetic fields. We neglect the acceleration of the electrons 
due to electric fields in our computations. A discussion of 
the effect of electric fields can be found in Hededal 2005, 
who conclude that the effect is that of flattening the low 
frequency slope. This effect will be important in the mod- 
eling of the prompt bursts but not in afterglows, where the 
resultant spectrum will be dominated by the limb. 

• We also assume that there is no correlation between the 
electron energy and the magnetic field locally. Such assump- 
tion is customary also for synchrotron computations. 

• Finally, the limit 8 « 1 is assumed 

We find that depending on the orientation of the line 
of sight with respect to the shock front the jitter spectrum 
is different. For shocks observed head on, the spectrum is 
peaked, with a steep low energy slope F(y) oc v. For shocks 
observed edge on, the spectrum at the left of the peak is flat 
down to the self absorption frequency or to the frequency at 
which scattering off the magnetic field inhomogeneities be- 
comes important (Fleishman 2006b). We applied these re- 
sults to standard GRB cases, for the prompt and for the 
afterglow emission. We find that, in addition to some of 
the observations already discussed in the previous literature 



8 Workman et al. 




Figure 1. Jitter emissivity vs. frequency at different angles between the bulk velocity vector and the line of sight. The angles plotted 
arc 6' = 0° (solid line), 5° (long-dashed line), 15° (dashed line), 30° (dot-dash line), 45° (3 dot-dash), and 90° (dotted line). All the 
angles are measured in the comoving frame. 



(Medvedev 2000, 2006), jitter radiation can explain X-ray 
excesses in the prompt spectrum and the lack of correlation 
between optical and high energy radiation. This is due to 
the fact that the optical radiation lies in the regime where 
the spectral slope depends on the orientation angle, while 
the high energy emission depends only on parameters that 
do not evolve with time. We also find that jitter radiation 
explains naturally the presence of spectral lags (Norris et al. 
2001). 

In the afterglow regime, the optical and X-ray spec- 
tra are very similar to those of synchrotron radiation. At 
frequencies below the peak, differences are instead clearly 
visible. First, the spectrum is flat {F(v) oc v° instead of 



F(y) oc v 1 ^ 3 ); second, the self absorption frequency is at 
higher frequencies (see also Medvedev, et al. 2007 who also 
show that self absorption frequency has a different tempo- 
ral evolution in jitter and synchrotron). Understanding the 
implications of those differences is not simple and a proper 
fit has to be performed (Morsony et al. in preparation). It is 
likely that jitter spectra fit to afterglow data will give differ- 
ent results in term of the properties of the ambient medium, 
since radio observations are always very important in con- 
straining the density of the interstellar material. A more 
detailed discussions will be possible only after a formal fit 
has been performed. 
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Figure 2. Absorption coefficient a vs. frequency at different angles between the bulk velocity vector and the line of sight. The angles 
plotted are 0' = 0° (solid line), 5° (long-dashed line), 15° (dashed line), 30° (dot-dash line), 45° (3 dot-dash), and 90° (dotted line). 
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Figure 5. Jitter radiation light curves for a single thin relativistic shell at 1 MeV (dot-dash line), 100 keV (dashed line), 10 keV (solid 
line) and 1 eV (dotted line). 
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APPENDIX A: DERIVATION OF (|W W ,| 2 ) 

The following derivation is directly derived from Medvedev 
(2006) and Fleishman (2006b) and is included for the sake 
of completeness. First, let us define a parameter 8 which de- 
fines the ratio of deflection (a) due to Lorentz forces and 
beaming (A#) due to relativistic effects experienced by a 
particle moving with Lorentz factor 7 in a small scale, ran- 
dom magnetic field with a typical correlation length, Ub and 
Larmor radius, p e . 

s = T 1 -tt; ( A1 ) 

For values of S <JC 1 the particle trajectory is nearly a straight 
line (with small perpendicular motions introduced by the 
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Figure 6. Simulated BASTE 3:2 channel hardness ratio vs. time. HR3/2 is defined as the ratio of the number of photons between 110 
and 320 keV to the number of photons between 55 and 110 keV. 



Lorentz forces). In this case the angle averaged spectral en- 
ergy, neglecting plasma dispersion, is given by Landau & 
Lifshitz (1971, section 77, p. 215): 



dW _ e 2 oj 

dljj 27TC 3 



/2 7 2 



1 - 



al'7 2 



dJ (A2) 



where w w / is the Fourier transform of the transverse accel- 
eration of the particle due to Lorentz forces (w = Ft /7m). 
Fourier transforming the acceleration field along a particles 
trajectory w(ro + vi, t) yields 



= (27r)~ 3 J wn M S(cu' - Q + k ■ v) e*' r ° dQdk. (A3) 

Squaring the above result and averaging over vol- 
ume using the results (|w^/| 2 } = V~ 1 J |w w / | 2 dro and 
/ e '< k - kl >' r ° dr = (27r) 3 <5(k - k x ) yields 

<|w^| 2 ) = {2n)-' i V- 1 J |wn :k | 2 <5(a/-n + k- v)d«dk.(A4) 

In the absence of electric fieldfl the Lorentz acceleration 
is given by (e/ymc)v x B which, in tensor notation, is 



w u / = (2tt)- 4 / e ! "' t ^(e- i(nt - kro - k vt) wn,kdf2dk) 



2 See § 2.1 for a discussion on the implications of this assumption. 
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Figure 7. The total afterglow power emitted at 1000 s for jitter (solid line) and synchrotron (dashed line) radiation mechanisms under 
identical conditions. See text for details on the computations. 



pression results in 



v-yBp). After simplification this ex- 



|wn, k | 2 = (ev/'yrnc) 2 (5 a f3 - v 2 v a v/3) Bq M B^. 



(A5) 



where Bq ^B^f k ^ s the Fourier Transform of the field corre- 
lation tensor 

B%, k B^ k = TVK olP (n,k) = TvJe t{nt - k r) K a p(r,t)drdt, (A6) 

where T can be interpreted as the period or duration of 
emission, V is the volume integrated over and K a p(r,t) is 
the second order correlation tensor of the magnetic field. 
Using the results above and assuming a time independent 
magnetic field, equations (|A4[I and (|A5[) reduce to 



'l 2 ) = (27rvy 



J |wk| 2 <5(w' + k • v) dk, 



(A7) 



|w k | 2 = (ev/^mcf(8 a p - v 2 v a vp) VTK a p(k). (A8) 

The form of K a j}(k.) used in this paper is taken from nu- 
merical simulations of the Weibel Instability and is given 
by 



K a p(k) = C(8 a /3 - n a np)f z (k\\)f X y{kA 



(A9) 



where n is the normal to the shock front, C is propor- 
tional to the mean square magnetic field {B 2 }, and 
and fxy(k±) describe the structure of the magnetic field 
parallel and perpendicular to the shocks normal. Inserting 
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the above into equations (|A4|) and (|A5[) and simplifying 

(S a /3~v a vp/v 2 )(S al 3~n a n/3) = l + (n a v a ) 2 /v 2 = 1 + cos 2 9', 
finally gives the result (for v ~ c; Medvedev 2006) 

<|w.f> = (-^-) 2 ^(l + cos 2 e')x 

x // z (fc||)/^(fcx) ( S(Lj' + k-v)rffc||rf 2 fcx. (A10) 



